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Evidence is reported for a narrow structure near the J/'tp(t> threshold in exclusive B^ J /ipcpK'^ 
decays produced in pp collisions at -/H — 1.96 TeV. A signal of 14 ± 5 events, with statistical 
significance in excess of 3.8 standard deviations, is observed in a data sample corresponding to an 
integrated luminosity of 2.7 fb~^, collected by the CDF II detector. The mass and natural width 
of the structure are measured to be 4143.0 ± 2.9(stat) ± 1.2(syst) MeV/c^ and 11.7t|:^(stat) ± 
3.7(syst) MeV/c^ 

PACS numbers; 14.40.Gx, 13.25.Gv, 12.39.Mk 
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Heavy quarkonium spectroscopy provides insight into 
strong interactions that are not precisely predictable by 
QCD theory. The recently discovered states that have 
charmonium-like decay modes [H |3, S, 0] but are dif- 
ficult to place in the overall charmonium system have 
introduced new challenges. The possible interpretations 
beyond quark-antiquark states (qq) such as hybrid (qqg) 
and four-quark states (qqqq) have revitalized interest in 
exotic mesons in the charm sector [H, @, 0, S)- An im- 
portant tool in unraveling the nature of the states in the 
charmonium-mass region is the exploration of states in 
diverse channels. First, the J/ipff' final state, with posi- 



tive C-parity, and two J 



PC _ 1 — 



vector mesons (VV), 



is a good channel for an exotic meson search. The dis- 
covery of the X{3872) (proposed as a four-quark state 
(ccqq) 0,11) and 1^(3930) [f, both decaying into VV 0, 
suggests searching for other possible VV states [l^ . Sec- 
ond, the observation of F(3930) near the J/^poj thresh- 
old motivates searches for similar phenomena near the 
J I'ipi' threshold. Third, the observation of the y(4260), 
a potential hybrid candidate Q, leads to an expecta- 
tion of a triplet of hybrid states J^'~^ — (0,1,2) ^ to 
lie nearby in mass [ill [l^ . which would be accessible 
in the J/'4'4' channel. Finally, other possibilities such 
as glueballs [l3| and nuclear-bound quarkonium [l3| also 
motivate a search in this channel. The J/ip(f> channel 
is accessible in the decay mode _B+ J/tfj(j)K~^, which 
has been observed |lj|. However, to date no results have 



been reported for substructure in the J /ipfj) channel. 

In this Letter, we report an investigation of the J/'4'4> 
system produced in exclusive — > J/ip(j)K^ decays 
with J/tp ^ /^^M^ and (p ^ K^K^ . The search in 
exclusive decays is more sensitive than an inclu- 
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sive search since the additional B~^ mass constraint on 
the J/il'4>K~^ system helps to reduce background. This 
analysis is based on a data sample of pp collisions at 
^/s = 1.96 TeV with an integrated luminosity of 2.7 fb~^ 
collected by the CDF II detector at the Tevatron. Charge 
conjugate modes are included implicitly in this Letter. 

The CDF II detector has been described in detail else- 
where [l5| . The important components for this analysis 
include the tracking, muon, and time-of-flight (TOF) sys- 
tems. The tracking system is composed of a silicon-strip 
vertex detector (SVX) surrounded by an open-cell drift 
chamber system called the central outer tracker (COT) 
located inside a solenoid with a 1.4 T magnetic field. 
The COT and SVX are used for the measurement of 
charged-particle trajectories and vertex positions. In ad- 
dition, the COT provides ionization energy loss infor- 
mation, dE/dx, used for particle identification (FID), 
while the TOF system provides complementary FID in- 
formation. The muon system is located radially out- 
side the electromagnetic and hadronic calorimeters and 
consists of two sets of drift chambers and scintillation 
counters. The central part of the muon system covers 
the pseudorapidity region I77I < 0.6 and detects muons 
with pt > 1.4 GeV/c [16|, and the second part cov- 
ers the region 0.6 < |r/| < 1.0 and detects muons with 
PT > 2.0 GeV/c. 

In this analysis, J ftp — s- /i^/^^ events are recorded us- 
ing a dedicated three-level dimuon trigger. The first trig- 
ger level requires two muon candidates with matching 
tracks in the COT and muon systems. The second level 
applies additional kinematic requirements to the muon 
pair candidate. The third level requires the invariant 
mass of the ji'^ pair to be within the range of 2.7 to 
4.0 GeV/c2. 

Offline reconstruction of i?+ J/tp(j)K^ candidates 
uses only tracks that pass standard CDF quality require- 
ments and which have been corrected for ionization en- 
ergy loss for the muon or kaon hypothesis, as appropri- 
ate. The -B+ J/ilj<j)K^ candidates are reconstructed 
by combining a J/tp —^ fJ-"^ candidate, a </) — *■ K~^K~ 
candidate, and an additional charged track. All five 
tracks must form a good quality 3D vertex, using a pri- 
ori requirements typical for B hadron reconstruction at 
CDF [18|. Preliminary event selection requires a, J/^p 
candidate reconstructed using opposite-sign muon candi- 
dates and a (p candidate formed from opposite-sign tracks 
to which we assign the kaon mass. Masses of vector me- 
son candidates must lie within ±50 MeV/c^ of the J/tp 
mass for muons or ±7 MeV/c^ of the (p mass for kaons. In 
the final B'^ reconstruction the J/ip is mass constrained, 
and the B'^ candidates must have pt > ^ GeV/c. 

To suppress combinatorial background, we use dE/dx 
and TOF information to identify all three kaons in the 
final state. The information is summarized in a log- 
likelihood ratio (LLR), which reflects how well a candi- 
date track can be positively identified as a kaon relative 
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FIG. 1: (a) The mass distribution of J/tpcpK'^; the solid line 
is a fit to the data with a Gaussian signal function and flat 
background function, (b) The sideband-subtracted mass 
distribution of K'^K~ without the mass window require- 
ment. The solid curve is a P-wave relativistic Breit-Wigner 
fit to the data. 



FIG. 2: (a) The Dalitz plot of m'^{(t)K+) versus m^{J/ip(j>) 
in the B'^ mass window. The boundary shows the kine- 
matic allowed region, (b) The mass difference, AM, between 
^~ K" and in the mass window. The dash- 

dotted curve is the background contribution and the red solid 
curve is the total unbinned fit. 



to other hadrons [19]. In addition, we require a mini- 
mum Lxy{B^) for the i?^ J/ip(j)K^ candidate, where 
Lxy{B^) is the projection onto pt{B^) of the vector con- 
necting the primary vertex to the i?"*" decay vertex. The 
primary vertex is determined for each event using prompt 
tracks. 

The Lxy{B^) and LLR requirements for i?+ 
J/ ip(l)K^ are then chosen to maximize S/y/S + B, where 
S is the number of B^ J/tlj(j>K^ signal events and 
B is the number of background events in the J /ijxpK^ 
mass range of 5.0 to 5.6 GeV/c^ in the data. The val- 
ues of S and B are determined from an unbinned log- 
likelihood fit to the mass spectrum of J/4>4>K'^, for a 
given set of values of Lxy{B^) and LLR. A Gaussian 
function is used to represent the 5+ J/ijj(j>K'^ signal, 
where the mean value of the Gaussian is fixed to the _B+ 
world-average mass value The B~^ mass resolution 
is fixed to the value 5.9 MeV/c^ obtained from Monte 
Carlo (MC) simulation A linear function is used to 
model the background in the fit. The requirements ob- 
tained by maximizing S / \^ S + B are L^y {B^ ) > 500 /im 
and LLR > 0.2. In order to study the efficiency of 
the Lxy{B~^) and LLR selections, we also reconstruct 
5+ — > J/tpK^ and J/ipff) as control channels. We 

select approximately 50 000 B+ J/^K+ and 3000 
— !■ J/ip4' events by applying similar requirements 
as for the J/ipcfiK'^ channel but without the Lxy{B^) 
and LLR requirements. The efficiency for PID with the 
LLR > 0.2 requirement is approximately 80% per kaon 
and is reasonably flat as a function of kaon pr; the ef- 
ficiency for Lxy{B^) > 500 /xm is approximately 60%, 
based on the B^ — > J/tpK^ control sample. 

The invariant mass of J/ip<j)K^ after the Lxy{B^) and 
LLR requirements and J/tp and (j) mass window require- 
ments is shown in Fig. [TJa) . A fit with a Gaussian signal 
function and a flat background function to the mass spec- 
trum of J/tp(l)K~^ returns a B~^ signal of 75 ± lO(stat) 
events. We select B~^ signal candidates with a mass 



within 3cr (17.7 MeV/c^) of the nominal B+ mass; the 
purity of the B"*" signal in that mass window is approxi- 
mately 80%. 

The combinatorial background under the B~^ peak 
includes B hadron decays such as i?° — > tp{2S)(l) 
J/ijjTT^Tr^(f), in which the pions are misidentified as kaons. 
However, background events with misidentified kaons 
cannot yield a Gaussian peak at the B'^ mass consistent 
with the 5.9 MeV/c^mass resolution. The kinematics are 
such that for the hypothesis B~^ J /ipK'^ K~ K'^ , only 
events with real kaons can produce the observed Gaus- 
sian signal. Thus, with the B"*" mass window selection 
the sample consists of real B^ J/tpK^K~K^ decays 
over a small combinatorial background. 

Figure [Hb) shows the invariant mass distribution of 
K^K~ pairs from jX^ ii~ K~ candidates within 
±3(7 of the nominal B"*" mass. The spectrum shown in 
this figure has had the sidebands subtracted, but the (j) 
mass window selection has not been applied. By fitting 
the K^K~ mass spectrum to a P-wave relativistic Breit- 
Wigner (BW) function 2l| convoluted with a Gaussian 
resolution function with the rms fixed to 1.3 MeV/c^ ob- 
tained from simulation, we obtain a mass of 1019.6 ± 0.3 
MeV/c^ and a width of 3.84 ± 0.65 MeV/c^with 
probability of 28%, consistent with the world-average 
values for the meson [l3|. The good fit indicates 
that after the ±7 MeV/c^ selection on the (j) mass win- 
dow, the B^ J/tfjK~^K~ final state is well de- 
scribed as J/ijj(j)K^, with negligible contributions from 
J/ipfo{980)K+ or J/^K+K~K+ phase space. 

We examine the effects of detector acceptance and se- 
lection requirements using B^ J/tp(f>K'^ MC events 
simulated by phase space distributions. The MC events 
are smoothly distributed in the Dalitz plot and in the 
J/tp(j> mass spectrum. Figure [2ja) shows the Dalitz plot 
of nn?{(f)K^) versus nn?{J /ipcf)), and Fig. [2l^b) shows the 
mass difference, AAf = m(fi'^ ^~ K~) — m{p,'^ij,~), 
for events in the B'^ mass window in our data sample. 
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We examine the enhancement in the AM spectrum just 
above J/^cj) threshold, using 73 events with AM < 1.56 
GeV/c^. We exclude the high mass part of the spectrum 
to avoid combinatorial backgrounds from misidentified 
5° ^ ip{2S)(j) (J/'07r+7r-)0 decays. 

We model the enhancement by an S'-wave relativistic 
BW function convoluted with a Gaussian resolution 
function with the rms fixed to 1.7 MeV/c^ obtained from 
MC, and use three-body phase space [l7| to describe 
the background shape. An unbinncd likelihood fit to the 
A Af distribution, as shown in Fig. [H^b) , returns a yield of 
14 ±5 events, a AM of 1046.3 ±2.9 MeV/c^ and a width 
of ll-ltll MeV/c^. We also fit the AM distribution to 
a single Gaussian with rms given by the mass resolution 
(1.7 MeV/c^), plus phase space background, to test the 
hypothesis that the structure has zero width. The sta- 
tistical significance for a non-zero width determined by 
the log-likelihood ratio between these two fits is 3.4a-, 
indicating a strong decay for this structure. 

We use the log-likelihood ratio of —2\r^{Co / Cmax) to 
determine the significance of the structure at the 
threshold, where Cq and Cmax are the likelihood val- 
ues for the null hypothesis fit and signal hypothesis fit. 
The —2\n{Co / Cmax) value is 5.3 for a pure three-body 
phase space background shape assumption. To estimate 
the probability that background fluctuations alone would 
give rise to signals as significant as that seen in the data, 
we simulate AM spectra based on the background dis- 
tribution alone, and search for the most significant fluc- 
tuation in each spectrum in the mass range of 1.02 to 
1.56 GeV/c^, with widths in the range of 1.7 (resolution) 
to 120 MeV/c^ (ten times of the observed width). From 
these spectra we obtain the distribution for the quantity 
—2\a.{CQ I Cmax) in pure background samples, and com- 
pare this with the signal in the data. We performed a 
total of 3.1 million simulations and found 29 trials with 
a —2\ti[Cq / Cmax) value greater than or equal to the value 
obtained in the data. The resulting p-value is 9.3 x 10"^, 
corresponding to a significance of 4.3(t. Thus, the signifi- 
cance is decreased from a simple estimate of 5.3(T to 4.3(7 
by taking into account the absence of a prior prediction 
for the mass and width [23] . 

In the analysis described above, we assumed that the 
backgrounds to the BW signal from both 5+ J/il}(f)K^ 
decays and combinatorial events in the mass window 
are described by three-body phase space. After making 
the AM < 1.56 GeV/c^ selection, we perform a fit to 
the J/ilj<j)K'^ mass spectrum with a Gaussian i?+ sig- 
nal and an empirical linear background shape; from this 
fit we determine that the 73 events within the final 
mass window include 15 ± 1 combinatorial background 
events by scaling it from the entire J/ip(j}K^ mass range. 
We model the i?+ events using three-body phase space 
as above, but use a flat spectrum to describe the combi- 
natorial events. This increases the average background 
level at small AM and reduces the yield by one event. 



The y/ —2\n{Co / Cmax) value with this modeling of back- 
ground is 4.8. We performed a total of 1.1 million simu- 
lations and found 99 trials with a —2hi{Co / Cmax) value 
greater than or equal to the value we obtained in the 
data. The p-value determined by this MG simulation is 
9.0 X 10^^, about 3.8cr significance. 

The mass of this structure is 4143.0 ± 2.9 MeV/c^ after 
including the world-average Jf-ip mass. To study the sys- 
tematic uncertainties of the mass and width, we repeat 
the fit to the AM distribution while varying the back- 
ground shapes as described above, and separately switch- 
ing to a non-relativistic BW for signal. The largest de- 
viation from the nominal values are 1.2 MeV/c^ for AM 
and 3.7 MeV/c^ for the width. Therefore we assign a sys- 
tematic uncertainty of 1.2 MeV/c^ to the mass and 3.7 
MeV/c^ to the width. 

There is a small cluster of events approximately one 
pion mass higher than the first structure, located around 
1.18 GeV/c^ in Fig. [IJb). However, the statistical sig- 
nificance of this cluster is less than 3a. To investigate 
possible reflections, we examine the Dalitz plot and pro- 
jections into (j)K~^ and J/ijjK'^ spectrum. We find no 
evidence for any other structure in the ^iC + and J / ipK^ 
spectrum; the only structure [i.e. ^2(1770)] that has 
been claimed in the (j)K~^ spectrum by previous experi- 
ments is too broad to alter our analysis ^231] . 

In summary, the large sample of B~^ J/ip(l)K~^ de- 
cays (75 events) enables us to search for structure in the 
J/ipcj) mass spectrum, and we find evidence for a narrow 
structure near the J/ipej) threshold with a significance in 
excess of 3.8cr. Assuming an S'-wave relativistic BW, the 
mass and width of this structure, including systematic 
uncertainties, are measured to be 4143.0 ± 2.9(stat) ± 
1.2(syst) MeV/c2 and 11.7t|^(stat) ± 3.7(syst) MeV/c^, 
respectively. It is well above the threshold for open charm 
decays, so a cc charmonium meson with this mass would 
be expected to decay into an open charm pair dominantly 
and to have a tiny branching fraction into J/tp(f> [5[. 
Thus, this structure does not fit conventional expecta- 
tions for a charmonium state. We note that this structure 
decays to J/ip4> just above the Jip(j) threshold, similar to 
the y(3930) [2], which decays to J/ipuj near the J/ifju) 
threshold. We therefore term it y(4140). 
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